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pH I We study the statistical and geometrical properties of the potential temperature 

^ ', (PT) field in the Surface Quasigeostrophic (SQG) system of equations. In addition to 

> : 

fSJ , extracting information in a global sense via tools such as the power spectrum, the g- 

O ; 

CD I beta spectrum and the structure functions we explore the local nature of the FT field 

I by means of the wavelet transform method. The primary indication is that an initially 

o : 

, smooth PT field becomes rough (within specified scales), though in a qualitatively 

I sparse fashion. Similarly, initially ID iso-PT contours (i.e., PT level sets) are seen 

I to acquire a fractal nature. Moreover, the dimensions of the iso-PT contours satisfy 

X ; 

^ , existing analytical bounds. The expectation that the roughness will manifest itself in 

the singular nature of the gradient fields is confirmed via the multifractal nature of the 
dissipation field. Following earlier work on the subject, the singular and oscillatory 
nature of the gradient field is investigated by examining the scaling of a probability 
measure and a sign singular measure respectively. A physically motivated derivation 
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of the relations between the variety of scahng exponents is presented, the aim being 
to bring out some of the underlying assumptions which seem to have gone unnoticed 
in previous presentations. Apart from concentrating on specific properties of the 
SQG system, a broader theme of the paper is a comparison of the diagnostic inertial 
range properties of the SQG system with both the 2D and 3D Euler equations. 

PACS number 47.27 



2 



1. Introduction 

In the Quasigeostrophic (QG) frameworkS, a simplification of tlie Navier Stokes equations for 
describing the motion of a stratified and rapidly rotating fluid in a 3D domain, there are two 
classes of problems that immediately come to attention. The first (Charney type) are the 
ones where attention is focussed on the interior of the domain; the temperature is uniform 
along the boundaries and they play no dynamical role in the evolution of the system. The 
other class (Eady type) of problems lead to the Surface Quasigeostrophic (SQG) equations. 
The potential vorticity in the 3D interior is forced to be zero and the dynamical problem is 
controlled by the evolution of the potential temperature at the 2D boundaries. Working with 
a single lower boundary (assuming all fields to be well behaved as z oo), the equations 
making up the SQG system can be expressed ad^,!, 

dtO + u'diO = Dte = ,z = 

e = d,iij = (-A)^ 

where V^?/^ = 0, 2; > and u = V^ip (1) 

Here 9 is the potential temperature (it is a dynamically active scalar due to the coupling of 
6 to Tp) ^ Tp \s the geostrophic streamfunction, u the geostrophic velocity, V"*- = {—dy, d^), A 
is the horizontal Laplacian, is the full 3D Laplacian, the operator (—A) 2 is definedS in 
Fourier space via {—A)'^ip{k) = \k\%l)(k) and i = x,y. Recall that the 2D Euler equations 
(for an incompressible fluid) in vorticity form are, 

dti + u'dii = Dt^ = 

^ = Atjj with u = V^i) (2) 

where ^ is the vertical component vorticity. The similarity in the evolution equations for 
6 and ^ has been explored in detail by Constantin, Majda and Tabak0,i. It can be seen 
that the structure of conserved quantities in both equations is exactly the same. To be 
precise, just as f{^)J are conserved by the 2D Euler equations similarly f{6)J ipO are 



conserved in the SQG system. The basic difference in the above two systems is the degree 
of locahty of the active scalar. For the 2D Euler equations the free space Green's function 
behaves as In (r) implying a ^ behaviour for the velocity field due to point vortex at the 
origin. In contrast, in the SQG equations the free space Green's function has the form ^ 
implying a much more rapidly decaying ^ velocity field due to a point "PT vortex" at the 
origin!. Or, in Fourier space one has ^ = \k\'^'ijj and 6 = \k\'il) for the 2D Euler and SQG 
equations respectivelyi. Hence the nature of interactions is much more local in the SQG 
case as compared to the 2D Euler equations. 

Studying the properties of active scalars with different degrees of locality would be an 
interesting question in its own right! but the specific interest in the SQG equations comes 
from an analogy with the 3D Euler equations. This can be seen by a comparison with the 
3D Euler equations, which in vorticity form read, 

[^Y = d^v^co^ (3) 

where w is a divergence free velocity field, uj{= V x -u) is the vorticity and i,j = x,y,z. 
Introducing V = a "vorticity" like quantity for the SQG system which satisfies (dif- 

ferentiating Eq. (|l]) and using incompressibilty) , 

l^y = d^uw^ (4) 

Identifying V in Eq. (^) with uj in Eq. it can be seenH that the level sets of 6 are geomet- 
rically analogous to the vortex lines for the 3D Euler equations. Similar to the question of 
a finite time singularity in the 3D Euler equations (which is thought to be physically linked 
to the stretching of vortex tubes), in the SQG system one can think of a scenario where the 
intense stretching (and bunching together) of level set lines during the evolution of a front 
leads to the development of shocks in finite time. The issue of treating the SQG system as a 
testing ground for finite time singularities has generated interest@,i,0,i,lHl in the mathematical 
community and the reader is referred to the aforementioned papers for details regarding this 



issue. 



In view of the similarities between the SQG and 2D Euler equations and the level set 
stretching analogy with the 3D Euler system, it is natural to inquire into the statisti- 
cal/geometrical properties of the SQG active scalar within an appropriately defined "inertial 
range". The broad aim is to compare these properties with the large body of work avail- 
able for the 2D and 3D Euler equations. In the second section we examine the PT field 
via global (power spectrum, structure functions, {P,g{P)) spectrum) and local (wavelets) 
methods. One of the few rigorous estimates that exist in fluid turbulence is that for the 
level set dimensions. The extraction of these dimensions and their agreement with analytical 
bounds is demonstrated. The third section is devoted to the examination of the dissipation 
field, the generalized dimensions of a measure based upon the dissipation field are calculated 
and commented upon. In the fourth section we focus our attention on the gradient fields, a 
simple derivation of the relation between the variety of scaling exponents is presented and 
the underlying assumptions are clearly stated. The failure of the cancellation exponent is 
demonstrated and a simple example is presented so as to put some of the ideas in perspec- 
tive. 

2. The Potential Temperature Field 

A. The Power Spectrum and the Structure Functions 

A pseudo-spectral technique was employed to solve Eq. (|l|) numerically on a 2048 x 2048 
grid. Linear terms are handled exactly using an integrating-factor method, and nonlinear 
terms are handled by a third-order Adams-Bashforth scheme (fully de-aliased by the 2/3 
rule method). The calculations were carried out for freely decaying turbulence. The initial 
conditions consisted of a large-scale random field, specifically a random-phase superposition 
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of sinusoids with total wavenumber approximately equal to 6, in units where the gravest 
mode has unit wavenumber. Potential temperature variance is dissipated at small scales 
by vV'^Q diffusion. Based on the typical velocity V and scale L of the initial condition, 
one may define a Peclet number UL/v. The calculations analyzed here were carried out 
for a Peclet number of 2500. After a short time, the spectrum develops a distinct inertial 
range. As time progresses, energy and 9 variance are dissipated at small scales, the am- 
phtude decreases, and the effective Peclet number also decreases. After sufficient time, the 
flow becomes diffusion-dominated and the inertial range is lost. Analysis of other cases, not 
presented here, indicates that the results are not sensitive to the time slice or the Peclet 
number, so long as the Peclet number is sufficiently large and the time slice is taken at a 
time when there is an extensive inertial range. 

The mean ID power spectra from different stages of evolution can be seen in Fig. |l|. 
As these are decaying simulations the structure in the PT field is slowly dying out. The 
resulting increase in smoothness of the PT field can be seen via the roll off of the spectrum 
during the later stages. In spite of this a fairly clean power law is visible for a sizeable 
"inertial range" (other runs with large scale initial conditions posessing various amounts of 
energy show similar behavior). We choose to concentrate on the particular stage which has 
the largest inertial range. The 2D power spectrum for this stage can be seen in Fig. § and 
a snapshot of the PT field itself can be seen in Fig. ^. Interestingly the 2D power spectrum 
seems to roll off at larger wavenumbers as compared to the ID spectrum. In this stage the 
spectral slope (from the ID spectra) between the scales 256 to 8 (the scales are in terms of 
grid size) is ~ —2.15 (the other runs also showed slopes steeper than -2). The slope from the 
2D spectrum is ~ —2.11 (due to the early roll off of the 2D spectrum, this slope is extracted 
between the scales 128 to 8). Previous decaying simulation^ obtained values near —2 and 
seem to be consistent with our observations. A slope as steep as this suggests that the field 
being examined is smoother than expectations from a similarity hypothesis (which yields a 
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^ slopei). A closer look indicates (Fig. |^) that the field is composed of a small number 
of "coherent structures" superposed upon a background which has a filamentary structure 
consisting of very fine scales. This immediately brings to mind the studies on vorticity 
in decaying 2D turbulence0,0,li wherein a similar coherent structure/background picture 
was found to exist. Further analysis indicated that the vorticity field possessed normal 
scaling whereas a measure based upon the gradient of the vorticity (precisely the enstrophy 
dissipation) was multifractal. To proceed in this direction we introduce the generalized 
structure functions of order q G 3?"*", 

Sg{r) =<\e{x + r)-e{x)\'^ > (5) 

Here < . > represents an ensemble average. The directional dependence is suppressed due 
to the assumed isotropy of the PT field. Scaling behaviour in the field implies that one can 
expect the generalized structure functions to behave aslll, 

S,ir) = C^{q)\A9{L,)n^f ; n < r < L, (6) 

where (g are the generalized scaling exponents, Ci{q) is of order unity for all q, |A^(Li)| is 
the absolute value of the difference in 6 over a scale Li. ri and Li are the inner and outer 
scale (8 and 256 respectively) over which the power law in the spectrum was observed. 



If the field being examined is smooth at a scale r then the gradient at this scale would 
be finite and as a consequence Cq = Q (due to the domination of the linear term in the Tay- 
lor expansion about the point of interest) ie. the scaling would be trivial. Conventionally 
normal scaling is a term reserved for linear (g and any nonlinearity in (g is referred to as 
anomalous scaling. In 2D turbulence the velocity field is known to be smooth for all time 
if the initial conditions are smooth0 and henceS Cq(velocity)= q. Also, as mentioned, from 
the analysis of the vorticity field0 the scaling exponents for the vorticity structure functions 
were found to depend on g in a linear fashion. Plots of log(5'g(r)) Vs log(r) for the PT can 
be seen in the upper panel of Fig. ^. In all cases the scaling is valid upto r ~ 128, using 

7 



these plots we extracted (g which are presented in the lower panel of Fig. ^. It is seen 
that the scaling is anomalous and in fact a best fit to the scaling exponents is of the form 
Cq = Aq^, q>0 with B = 0.82. 

For the special case of g = 2, one can in principle relate the scaling exponent (2 to 
the slope of the power spectrum (n) via0 n = — (1 + (^2)- This relation is only valid for 
—3 < n < —1 (note that this does not prevent the spectral slope from being steeper than 
—3; it just implies that (2 saturates at 2 for smooth fields and the particular relation between 
(2 and n breaks down). In our case (2 = 1-05 so the predicted spectral slope is n = —2.05 
which is near the observed mean value of —2.15 (or —2.11 from the 2D spectrum). Even 
though the scaling exponents give an idea of the roughness in the field (anomalous scaling 
implying differing degrees of roughness) there is a certain unsatisfactory aspect about the 
structure functions, namely, there is no estimate of "how much" of the field is rough. The 
following subsection aims to address this very issue. 



B. The {I3,g{l3)) spectrum 

In scaling literature the roughness of a field is specified by means of an exponent (3 (> 0) 
defined a s§, 

|0(a; + r) -0(s)| ~ |r|^(^) (7) 

Here /3 is a function of position and it refers to the fact that the derivative of 9 will be 
unbounded as r — if /? < 1. As mentioned previously there is a lower scale associated 
with the problem so technically nothing is blowing up and in effect (3 < 1 represents the 
regions where the derivative will be large as compared to the rest of the field. Note that 
9 itself cannot be singular due its conserved nature. The focus is on whether an initially 
smooth 9 field becomes rough so as to cause the gradient fields to experience a singularity. 
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It is clear that Eq. is by itself of not much use in characterizing 9 (as (3 depends upon 
position), in fact a global view of the specific degrees of rougness of 9 can be attained via 
the {P,g{P)) spectrum which we introduce next. 



An iso-/5 set is defined as the set of all x's where (3{x) = (3 and g{(3) is the dimension 
(to be precise g{l3) should be viewed in an probabilistic fashion^) of an iso-/3 set. The 
dimensions of iso-/5 sets are derived by Frischlli. Briefly, at a scale r the probability of 
encountering a particular value of (3 is proportional to r'^~3{P'> (where d = 2 for the 2D PT 
field and = 1 for ID cuts of the PT field). By using a steepest descent argument in the 
integral for the expectation value of \9{x + r) — 9{x)\'^ one obtainsi,i, 

Cq = mm[q(3 + d- g{(3)] or g{(3) = max[g/5 + d-Q (8) 

Hence given for a fixed q = using the first part of Eq. (^ , 

q. = ^ (9) 

^ df3 ^ ^ 

Denoting the value of P for which Eq. (H) is satisfied by /9* we have, 

g{P,) = d + q,p,-Cg, (10) 

Notice that, as the structure functions involve moments with positive q, they only pick out 
/3's such that /? < 1 and g{P) < d. The (/?, g{P)) spectrum seen in Fig. ^ hints at a hierarchy 
in roughness of the PT field. From the calculations we see that (3 G [0.26,0.6]. In Fig. ^ 
along with (g we have plotted the lines corresponding to (g = qPmin and (g = qPmax- As is 
expected these lines straddle the actual scaling exponents. For smaller values of q the scaling 
exponents are close to qPmax as Pmax has the largest associated dimension whereas for higher 
values of q the scaling exponents reflect the roughest regions and hence tend towards qPmin- 
Note that from Eq. (||) and Eq. ([lO|) the approximation (g = Aq^,B < l,q > implies 
that P becomes large and g{l3) — > d as g ^ 0. In essence the picture that emerges is that 
even though 9 appears to become rough (with differing degress of roughness), in fact, it is 
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the smooth regions that occupy most of the space. 



C. The roughness of ^ : a qualitative local view 

All of the previous tools, the power spectrum, the {[3,g{[3)) spectrum and the structure 
functions extracted information in a global sense. To get a picture of the actual positions 
where the signal may have unbounded derivatives, and to get a qualitative feel of the spare- 
seness of these regions, one has to determine the local behaviour of the signal in question. 
Recently the use of wavelets has allowed the identification of local Holder exponents in a 
variety of signals. The Holder exponents are extracted by a technique known as the wavelet 
transform modulus maxima (WTMM) method§,i°IJil,@. The modulus maxima refers to the 
spatial distribution of the local maxima (of the modulus) of the wavelet transform. In a 
crude sense the previous methods used ensemble averages of the moments of differences in 
9{x) as "mathematical microscopes" whereas in the wavelet method it is the scale parameter 
of the wavelet transform that performs this task. 



By using wavelets whose higher moments vanish one can detect sigularities in the higher 
order derivatives of the signal being analyzedil. Our previous results indicate that the 
PT field becomes rough, i.e., the first derivatives of the PT field should be unbounded. 
Hence, the particular wavelet we use is the first derivative of a Gaussian whose first moment 
vanishesEi,i.e., it picks up points where the signal becomes rough. The wavelet transform of 
a ID cut of the PT field can be seen in Fig. ^. The cone like features imply the presence 
of a rough spotHl. The modulus maxima lines are extracted from this transformed field and 
can be seen in Fig. ^. The value of the local Holder exponent can be extracted via a log 
plot of the magnitude of a particular maxima lineil. In essence, the presence of the cones in 
the wavelet transform indicate roughness in the PT field and the WTMM lines locate the 
positions of the rough spots. However, a closer look (the lower panel of Fig. |^) suggests, 
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qualitatively, that the rough regions are sparsely distributed (for eg. comparing with Fig. 8 
in Arneodo et. al0). This goes along with the observation in the last subsection that the 
rough regions were non space filling. In contrast, it has been found§,i,§ that the local 
Holder exponents {h{x)) associated with a ID cut of the velocity field in 3D turbulence 
satisfy —0.3 < h{x) < 0.7 for almost all x (the peak of the histogram being at |) implying 
that the velocity field in 3D turbulence is very rough or has unbounded first derivatives at 
almost all points. 



D. The Dimension of Level Sets 

Apart from being physically interesting, the level set (iso-^ contour) dimensions provide 
another link to the roughness of the scalar field. As the initial condition was a smooth 2D 
field, initially, any given level set [Eg^ for 6 = 6q) is a non space filling curve, i.e., the level 
set dimension {D{Ee^)) is one. In the case of 3D turbulence there exist analytical estimates 
of the scalar level set dimensions^. These estimates are seen to be numerically satisfied by 
a host of fields (both passive and active) in isotropic 3D and Magnetohydrodynamic (MHD) 
turbulenceii. It must be emphasized that these estimates come directly from the equations 
of evolution and are much more powerful than the phenomenological ideas we have been 
working with so far. The actual calculationil,!! is of the area of an isosurface contained in 
a ball of specified size, given a Holder condition on the velocity field, 

\u{x + y,t)-u{x,t)\ < ULijf (11) 
this area estimate leads to the bound, 

D{Ee) <2.5 + ^ (12) 

The bound in Eq. (|1^) is expected to be saturated!! above a certain cutoff scale. A valid 
extrapolationii for the level sets of PT in the SQG system reads, D^Eq) < 1.5 + (Ci)/2, 
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where Ci is given by Eq. (H) due to the equahty of the scahng exponents for the velocity 
field and the PT in the SQG system. 

In passing we mention that the analytical estimates are for the Hausdorff dimension 
whereas practically we compute the box counting dimension. We performed calculations for 
a variety of nominal (i.e., near the mean) level sets and Fig. |^ shows the log-log plots used 
in the calculation of the box counting dimension. As is the case for 3Di,i there appears to 
be a crossover in D{Eg), we find that the even though the dimension undergoes a change, 
the fractal nature seems to persist at smaller scales. For small r, D{Eg) 1.3 whereas for 
larger values of r we find D{E0) ^ 1.74 (see the caption of Fig. ^for details on the actual 
values of r). /^From our previous calculations Ci = 0.57, hence the analytical prediction is 
D{Eq) < 1.785. Furthermore, as the bound is expected to be saturated above the cutoff we 
see that the computed value of D^Ee) for large r is quite close to the analytical prediction. 
In all, apart from satisfying the analytically prescribed bounds (and indirectly indicating 
roughness in the PT field), the level set dimensions indicate that initially non space filling 
level sets acquire a fractal nature in finite time. 

3. Gradient Field Characterization 

The effect of the inferred roughness in the PT field will, as mentioned previously, be reflected 
in the singular nature of the gradient fields. In this section we proceed to examine a variety 
of fields which are functions of the PT gradient. The aim is to see if we can actually detect 
the expected singularities, and if so, to characterize them. 
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A. The Dissipation Field 



A physically interesting function of the gradient field is the PT dissipation, as it is connected 
to the variance of the PT field. The equation for the dissipation of 9 can be obtained by 
multiplying Eq. (|l]) by 9 and averaging over the whole domain, 

= -2u < {y9f > (13) 

Here i/(V^)^ is the dissipation field. Now consider the quantity fj,{x,r), 

^(x,r) = 4 / i^(V9?d'^x (14) 

r'^ JB{x,r) 

Physically this is the average dissipation in a ball of size r centered at x. Due to the 
smoothing via integration it is expected that fi{x, r) will be fairly well behaved through 
most of the domain with intermittent bursts of high values concentrated in the regions 
where the PT field is rough. The multifractal formalismii,@ ,§ provides a convenient way 
to characterize such "erratic" or singular measures. The techniquelHi consists of constructing 
a measure {fi with suitable normalization) and using its moments to focus on the singularities 
of the measure. The domain in which the field is defined is partitioned into disjoint boxes 
of size r and it is postulated (see for eg.0) that moments of /j, will scale as, 

^/i(x,r)'? ~ r"'-'?'^ (15) 

where the sum goes over all the boxes into which the domain was partitioned. Consider the 
set of points where the measure scales and denote the dimension of this set by /(a) (again 
a probabilistic view is more precise). By similar considerations as for the g{l3) spectrum it 
can be seen thatJli, 

Tq = min(gQ; — /(a)) and /(a) = max(ga — Tg) (16) 

The function Tg is further related to the generalized dimensionsS via Dg = Tg/[q — 1). 
Practicallyil a log-log plot of the ensemble average of ^{x,rY for different values of r gives 
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{Dg — d){q — 1). By knowing (g, Dg) one can use Eq. (|l^) to obtain the (a, /(a)) spectrum. 



Again, as we expect the scahng of any physical quantity to be restricted to a range of 
length scales (say to r^), it is preferable to work in terms of ratios of ^ where Vf, is the 
outer scale, is the inner scale and rt > r > Tq. In these terms the generalized dimensions 
can be expressed aJ^l, 



where /i(rb) is the measure on the outer scale and 6*2(5') is order unity for all q (a similar 
relation would hold if we replaced the {^{rbY) in the right hand side by < yu(r) but with 
different C2(g)'s). For r < Va the dissipation field is assumed to have become smooth via the 
action of viscosity. The /(a) spectrum and the generalized dimensions for the dissipation 
field can be seen in Fig. ^ The nontrivial behaviour of the generalized dimensions demon- 
strates that the dissipation field is singular (with different singularity strengths) within these 
scales. Also, from the calculations we find that /(I) < 2 which acts as a check that the 
dissipation in a finite volume is bounded (as /(I) is the dimension of the support of the 
singular regions)^. 

B. General Considerations 

The gradient squared nature of the dissipation field implied that it was positive definite. 
In principle one could conceive of fields which possess singularities but aren't sign definite. 
In order to characterize this possible sign indefiniteness Ott and co-workers studiediiJH 
sign singular measures and a related family of exponents called the cancellation exponents 
(Kq). It is immediately clear that singular nature is a prerequisite for the phenomenon of 
cancellation, the reason being that a nonsingular field is bounded and we can always add a 
constant so as to make the field positive definite and hence eliminate cancellation. Consider 




r (D,-d)(q-l) 



(17) 
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a sign indefinite ID field 6{x) (which will later be interpreted as a ID cut of the PT field) 
and construct, 

^\x,r) = {-j^ I — Ida:'); (18) 

V{x,r) = 1^ £^\^)dx'\ (19) 

As fi' is defined for the magnitude of the gradient field we can postulate (similar to Eq. (|1^ 
for yu), 

r (D'-d)(q-l) 

</z'(rr >=C3(g)(/(r,)^)(-) (20) 

' d 

Here d = 1 and D'^ are the generalized dimensions associated with fi'. As usual the scaling 
is restricted to a range of scales (rc < r < r^) and Vc is the scale below which fi' appears 
smooth. Formally, the entity r] after suitable normalization is a sign singular measure (see 
Ott et. al.ii for a rigorous definition). It was conjecturedil that the sign singular entity rj 
might also possess scaling properties in analogy with //', implying, 

< r]{ry >= C4{q){r)-< ] r > r^^ (21) 
Using Eq. (^0]) we can express Eq. ( ^1]) as, 

< 'q{ry >~< yu'(r)^ > (r)""' ; Vd > r > rcc > rc (22) 

Where + {D'^ — d){q — 1) and r^c is the lower oscillatory scale below which the 

derivative does not oscillate. We prefer to call Kg from Eq. (|22| ) the cancellation exponents 
as they directly reflect the difference in the scaling properties of ^ and | ^ | . 

Apriori there is no justification in assuming that the scale at which cancellation ceases 
(tcc) is the same as the scale at which fi' becomes smooth (rc). Consider the example where 
the derivative is a discrete signal composed of a train of delta functions (zero elsewhere) 
where the minimum separation between the delta functions is /. Furthermore let us assign 
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the sign of the delta functions in a random fashion. In this case the r^c = I, in fact, ki = 0.5 
due to the random distribution of the signs. But as the deha functions are supported at 
points we have Vc 0. As an aside we point out that if there is a maximum scale of sepa- 
ration between the delta functions (say L) then at scales greater than L we will see D'^ = d 
for fi'. The reason for pointing this out is to give a feel for fields that exhibit scaling, on 
one hand smooth fields have D'^ = d whereas on the other extreme random fields with small 
correlations also have D'^ = d (at scales larger than their correlation lengths). Fields with 
nontrivial scaling over a significant range are in effect random but with large correlation 
lengths. The reader is referred to Marshak et. al.il for a detailed examination of multi- 
plicative processes and the resulting characterization by structure functions and generalized 
dimensions. 



In order to get a unified picture of scaling in both the gradient and the field itself there 
have been attempts (for eg.Ei,0) to link Q, and D'^ to each other. The view that seems 
to have emerged is that there exist simple relations linking the various exponents and that 
these relations are valid under very general conditions. We present an alternate derivation 
of some of these relations which makes the implicit assumptions explicit. Proceeding from 
Eq. (p2D, assuming that ^ has integrable singularities we obtain, 

< \e{x + r) - e{x)\'' >=< 12' {ry > r'^-^^ (23) 
Substituting from Eq. ( pO]) and on comparing with Eq. (||) we get, 

Cg = iD'g-d)iq-l)+q-f,g (24) 
Writing Eq. (p3D for g = 1 we have, 

< \e{x + r) - e{x)\ >=< /i'(r) > r^^"' (25) 

which yields + Ci = 1- Now if we make a strong assumption regarding the uniform nature 
of the cancellation, it is possible to claim that Eq. ( ^Sf ) holds not only in average but on 
every interval, ie, 
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\e{x + r) - e{x)\ = fi'iry-''' (26) 
Raising this to the g*^ power and performing an ensemble average yields, 
< \e{x + r) - Oix)]" >=< fi'iry > r5(^-"i) (27) 
which implies, 

C, = (D;-rf)(g-l) + g(l-«:i) (28) 

The implications of Eq. (pSD are quite severe in that it shows to be dependent on D'^ and 
the knowledge of only the first cancellation exponent allows the derivation of one from the 
other. In general the scaling exponents of 6 and the generalized dimensions of fi' provide 
exclusive information. It is only in the presence of integrable singularities that one can link 
the two via Eq. (p^, furthermore the stronger relation (Eq. (^)) is valid under the added 
assumption of uniform cancellation. 



C. The Gradient of the PT and its absolute value 

We proceed to check if scaling is observed (as postulated) for the PT field gradient and its 
absolute value and whether one can extract the aforementioned exponents. In the upper 
panel of Fig. [1^ we show the log- log plots of < ^'{rY > Vs. r for different q. The scaling 
relations certainly appear to hold true (which was expected as they held for the dissipation 



field) . The generalized dimensions for /i' can be seen the lower panel of Fig. ITOl On the other 



hand the scaling for < ■r](r) > seen in Fig. |Tl| fails to exhibit a power law in r. Hence there 



is no meaningful way of extracting the cancellation exponents as they have been defined. 
Unfortunately this implies that the relations derived in the previous section (Eqs. (|2^) and 



(H)) cannot be used in this situation. In spite of this, we can see that as r decreases the 
tangent to log(< ?7(r) >) has a smaller slope which is consitent with the existence of a 
oscillatory cutoff at small scales. 
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4. Conclusion and Discussion 

In summary, we have found that the PT field in the SQG equations appears to become rough 
within a specified range of scales. Moreover, not only is there a heirarchy in the degree of 
roughness, the roughness is distributed sparsely in a qualitative sense. These conclusions are 
based on a combination of factors, namely, the algebraic power spectrum, anomalous scaling 
in the structure functions, a nontrivial {(3,g{(3)) spectrum, the nature of the WTMM map 
and the wrinkling of the PT level sets. The roughness in the PT field is expected to have 
an adverse effect on functions of the gradient field. This expectation is bourne out in the 
multifractal nature of the dissipation field. Also, the singular nature of the gradient field in 
combination with its sign indefinetness led us to examine a sign singular measure based upon 
the gradient field. The failure to observe scaling in the sign singular measure serves, in our 
opinion, as a reminder that most scaling arguments are postulated at a phenomenological 
level and the underlying basis of why scaling is observed in the first place is a fairly subtle 
and unsettled issue. Similarly the simple derivation of the relation between the variety of 
scaling exponents makes explicit some of the assumptions that are required for the validity 
of similar relations proposed in earher studies. 

Regarding the more general question we posed in the very beginning of this paper, i.e., 
where does the SQG active scalar stand with respect to both 2D and 3D turbulence, we have 
the following comments. With respect to the vorticity in the 2D Euler equations the cor- 
responding quantity in the SQG system is the PT. The coherent structure/fine background 
nature of the vorticity field0,0,i carries over qualitatively to the PT field. The vorticity 
structure functions which showed normal scaling in the 2D Euler system0 cross over to 
anomalous scaling in the SQG system. Our conjecture is that the stronger local interac- 
tions in the SQG system are responsible for this anomalous scaling. In both these systems 
the vorticity and PT respectively are conserved quantities and hence any singularities one 
might experience are actually in the gradient fields, as is seen in the multifractal nature of 
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the enstrophy dissipation in the Euler systemE3 and the PT dissipation in the SQG system. 



In the 3D Euler case, the PT from the SQG equations is analogous to the velocity field 
and from SQG corresponds to the vorticity of the 3D Euler equations. The roughness of 
the PT field is similar to the postulated roughness of the velocity field in the inertial range, 
a subtle difference being that the roughness in the PT appears to be sparse whereas indi- 
cations are that the roughness in the 3D velocity field is present almost everywhere@,0,il. 
We re-emphasize that in both these cases the roughness is restricted to a range of scales, 
i.e., no claim is made for an actual singularity in the corresponding gradients. Similarly 
the anomalous scaling of the PT follows that of the velocity field in 3D but again it is not 
as strong as in the 3D case. The general theory developedH for the deformation of scalar 
level sets in the 3D case is seen to carry over to the SQG equations. In essence the SQG 
equations follow the 3D Euler equations but in a somewhat weaker sense. This "weakness" 
is clearly manifested in the behaviour of the gradients. In the 3D Euler equations a sign 
singular measure constructed from the vorticity field shows good scaling properties and a 
cancellation exponent can be meaningfully extractedii, whereas in the SQG equations a 
similarly constructed entity lacks scaling. 



The results presented in this paper have been for the most part diagnostic, in that they 
characterize the nature of the roughness of the potential temperature field in SQG dynamics. 
Although we have exhibited anomalous scaling of the potential temperature fiuctuations, we 
do not have a theory accounting for the observed form of Q. Arriving at such a theory will 
be a major challenge for future work. Our results point efforts in the direction of considering 
the dichotomy between smooth fields within large organized vortices, and a rather sparse set 
at the boundaries of and between vortices which exhibits a greater degree of roughness. The 
diagnostic results also suggest a means for distinguishing between SQG and Euler dynamics 
in Nature, in cases where only a tracer field can be observed, as in the gas giant planets 
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(notably Saturn, Jupiter and Neptune, which exibit a rich variety of turbulent patterns). 
SQG dynamics should yield anomalous scaling corresponding to sparse roughness, whereas 
Euler dynamics should yield normal scaling. 
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FIGURES 
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Fig. 1. Power spectrum of the FT field for a variety of stages. The dashed hne (extracted from 



the stage which has the largest inertial range) has a slope of —2.15. 
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Fig. 5. The (/?, spectrum 
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Fig. 6. Upper Panel : A random ID cut of the PT field. Lower Panel : Its wavelet transform (the 
analysing wavelet is the first derivative of a Gaussian) , the cone like features indicate roughness in 
the analyzed signal. 
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Fig. 7. Upper Panel : The wavelet transform modulus maxima lines. The lower panel is a zoom 
into a particular section of the upper panel, this qualitatively indicates the sparseness of the rough 
spots. A log plot of the magnitude of the maxima line yields the local Holder exponent. 
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Fig. 8. log(A^(r)) Vs. log(r) for 3 different level sets. The dashed line is the best fit for 
4 < r < 60, the solid line is the best fit for 64 < r < 200. The box counting dimensions are 
(1.32, 1.8), (1.33, 1.77)and(l. 27, 1.64) for the small and large r regions for the 3 level sets respec- 
tively. 
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Fig. 9. The /(a) spectrum and Dg for the dissipation field 
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Fig. 11. log(< ry(r) >) Vs. log(r) 
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